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Abstract 

The recurrent relations between the eigenfunctions for GL{N,M) and GL{N — 1,M) 
quantum Toda chains is derived. As a corollary, the Mellin-Barnes integral represen- 
tation for the eigenfunctions of a quantum open Toda chain is constructed for the 
A -particle case. 
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1 Introduction 



Recently a new method of construction of the eigenfunctions for the periodic Toda chain have 
been introduced [IH . It rehes on generahzed Fourier transform expansion over eigenfunctions 
of an open Toda chain which coincide with the Whittaker functions of the GL(A^ — 1,]R) 
group i-i. 

In the original papers @-@] the Whittaker functions were constructed by purely algebraic 
methods in terms of the Iwasawa decomposition for the corresponding group. But it turns 
out that there is an alternative way to construct these functions directly on the level of 
i? -matrix formalism using essentially the integrable properties of the model. 

This paper was inspired by one remark of E. Sklyanin, that integral formula for the eigenfunc- 
tion of periodic Toda chain can be used to obtain some recurrent relations for -particles 
eigenfunctions of open Toda chain through the — 1 particles ones. Sklyanin's motivations 
based on a possibility to introduce a formal parameter to the Baxter equation |^ and then, 
putting it equal to zero, to transform the second order difference equation to the first order 
one. In our paper we find rigorous proof of this observation which is free of the limiting pro- 
cedure and find an universal way to construct "auxiliary" eigenfunctions staying completely 
in the framework of the i? -matrix formalism without any reference to an algebraic scheme 
developed in [^-lll- To compare with the well known results on eigenfunctions of open Toda 
chain ||^, |^ , our approach presents not only a new integral representation for the eigenfunc- 
tions but, together with generalized Fourier transformation, it gives a self-consistent method 
to solve the spectral problem for periodic Toda chain. We hope that this method can be 
applied for some other classes of integrable systems. 

2 The model 

We start with the i? -matrix formalism for the quantum periodic Toda chain 0. Let 

f X-Pn e"^" \ 
L„(A) = ^ (2.1) 

be the corresponding Lax operator where [2;„,pm] = ihSnm ■ The -particle monodromy 
matrix 

/ ^iv(A) B^{X) \ 

T^X) L,(A)...Li(A) ^ (2.2) 

V C^{X) D^{X) J 

satisfies the standard RTT relations with the rational -matrix. In particular, 

(A - ^ + ih)AMCN{X) = (A - fi)CMAM + ihA^WCM (2.3) 



The eigenfunctions for the periodic spectral problem have been constructed with the help 
of Weyl invariant function (xi, . . . , Xjv-i) = i^-yix) which is fast decreasing in the 
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regions Xk ^ 2:^+1, (A; = 1, . . . , — 1) and satisfies to equations 



iV-l 



C^{X)^^ = - e^- n (A-7^) i;^ (2.4) 



m=l 



Mij)^y = r^e-^^v^-^^e, (j = 1, . . . , iv - 1) 

(2.5) 

where is j -th basis vector in ^ . It is easy to see that the equation ( [^.5D is compatible 
with (2.41) due to commutation relation 



Actually, the function i/j-y{x) is an appropriate solution for A^ — 1 particle open Toda chain. 
Indeed, the operator Aj^_i{X) (arising from the A^— 1 particle problem) is the generating 
function for the Hamiltonians of the GL{N — 1, M) Toda chain. Using the obvious relations 
between the elements of the monodromy matrices Tf^{X) and Tjv-i(A) : 

A^{X) = (A-p,)A,_i(A) +e-^-C^_i(A) 

(2.6) 

C^(A) = -e^-A^_i(A) 
the equations ( |2.4| ) and ( |2.5|) can be written in the equivalent form 

N-l 

A^_i(A)^^= l[{X—fm)i^y (2.7) 

m=l 

CN-l{'Jj)i>'y = i'^lp-y-ihej {j = 1, . . . ,N - 1) 

(2.8) 

These equations fix (up to z/i -periodic common factor) the Weyl invariant Whittaker func- 
tion for G'L(A^— 1,R) group. In we choose the factor in such a way that i/j-^ is an entire 
function in 7 and the following asymptotics hold: 

i^f ~ l7,|— expj-— (A^-2)|7,|| (2.9) 

as |Re7j| 00 in the finite strip of complex plane. 

Remark 2.1 The function Cjv~i(A)'i/'-y is a polynomial in X of order N — 2 . Therefore, 
this polynomial is restored by their N — 1 values at given points 71, . . . ,'Jn-i ■ Hence, one 
obtains the interpolation formula 

A- 

c^_i(A)V'^ = Yl n -^r^ (2-10) 

Let us introduce the key object - the auxiliary function 

,T, / N dcf nV^-T. 7mjxjv (2.11) 
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where e is an arbitrary parameter. From ( p.7|) , (|2.1CI| ), and (|2.6| ) it is readily seen that this 
function satisfies to equations 

N-l N-1 N-1 , 

A,(A)vi>^,, = (-^ - e + 5^ 7m) n (^-^^■) ^-/.^ + ''"^ E n (2-12) 

C^(A)^^,, = - e^-- J] (A - 7,) (2.13) 

3 The problem 

Let the Weyl invariant Whittaker function '07i,...,7jv-i(^i) ■ ■ • ; ^n-i 

) for GL{N -1,R) Toda 

chain is given. The problem is to find the corresponding solution for GL{N, M) Toda 
chain using the above information, i.e. to construct the Weyl invariant Whittaker function 
4'\i,...,\Mi^ii • • • 5 ^n) satisfying to equations 

N 

^]v(A)^Ai,...,Aiv = n*^^ ~ ^'^^ ^Ai,...,Aiv (3-la) 
fc=l 

Civ(A„)?/^Al....,A^^ = '0Al,...,A„-^;^....,A^^ {u = 1, . . . , N) 

(3.1b) 

and obeying the asymptotics 

V-A^.^.A^ ~ |A„|— expj-— (iV-l)|A„|| (3.2) 

as I Re Art I — oo , in terms of the function i/j-y{x) . It is clear that the action of the operators 
An{X) and Cm{X) is nicely defined when acting on the auxiliary functions "^-^^^ixi, . . . , Xpf) . 
Therefore, it is reasonable to assume that the solution for GL{N, M) Toda chain is described 
by appropriate (generalized) Fourier transformation of the function \E'-y ,;(xi, . . . ,x^) . This 
is in complete analogy with the corresponding construction for the periodic case 



4 Main statements 



Theorem 4.1 Let \I'-y^(a;, Xjv) be the auxiliary function for N -periodic Toda chain, i.e. it 
is defined in terms of the Weyl invariant Whittaker function for GL(A^ — 1,]R) Toda chain 



according to Let A = (Ai, . . . , X^) G C be the set of indeterminates. Let 



^^{^) = {2nhr-\N-l)\X{ r(^) 



j<k 

N-l N 



(5(7i,...,7jv_i|Ai,...,Ajv) = Yl Yih-'^irr-T 

j=i k=i 



ih 



(4.1) 



(4.2) 
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Then the Weyl invariant Whittaker function for GL^N-M.) Toda chain is given by recurrent 
formula 

c 

where the integration is performed along the horizontal lines with Imjj > max^ {Im A^} . 

5 Proof of the Theorem 

First of all, the integral (^4.3|) is correctly defined. Indeed, the function 

N , 

g(7|Ai, . . . , A^) = n ^(^) ^^-^^ 



k=i 

obeys the asymptotics 



g(7;A)~|7r^/^exp{-^|7|} (5.2) 



as |Re7| — > 00 in finite horizontal strip while 



/i"'(7)~|7.r-'exp{^(iV-2)|7,|} (5.3) 



Hence, the integral in (|47^ ) is absolutely convergent due to asymptotics ( P^ 
Let us verify that the relation ( ^.laj ) holds. Using ( [2. 121 ) one finds 



^^(A)^Ai,.. •i-\jv (''"1 ) • • • 5 -^n) 
N N-l N~l 



/ (A-X]Afc + ^7„) JJ(A-7j)^ \j)Q{j;X)<iJ~f,Xi+...+XM{xi,---,XN)dj + 

^ k=l m=l j=l (5.4) 

N-l „ ^_ 

+ / n _ " /^'^(7)Q(7; A)^^-ite,,Ai+...+Ajv(3^i,---,3:^) ^7 

: 1 J 7.7 7m 



We shift the appropriate integrations 7-, + ih and use the functional equation 

7j - Im + ih 



+ ^h6,) = (-1) V^(7) n (5.5) 



The second integrand in ( |5.4| ) has no poles in any finite horizontal strip in the upper half- 
plane Im7j > maxfc{ImAfc} (all possible poles are cancelled by appropriate zeros of the 
function ). Moreover, the intergand is fast decreasing in this strip as Re7j — > ±00 . 
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As a consequence, the integral over the strip vanishes. Therefore, it is possible to deform 
the shifted contour to the original one. Hence, one arrives at the relation 



N N-l N~l 



^iv(A)l/'Ai,...,A 



AT 

c 

N-l 



I k=l m=l j=l 



(5.6) 



The function Q{-f; A) satisfies to equations 

N 

Uilj - h) Qir, A) = 2^g(7 + ^^e,; A) (5.7) 

k=l 

for any j = 1, . . . , N — 1 . Therefore, the relation (|5.6|) acquires the form 



c 

N-l N 



/( N N-l N-l 

^ I k=l m=l j=l 



(5.8) 



]=l k=l "i^j 



TV 

But the expression in curly brackets is nothing but the polynomial 11 (''^ ^ ^k) ■ Indeed, any 

k=l 

polynomial with the leading terms F{X) = + fiX^~^ + ■ ■ ■ can be uniquely restored by 
its values at any A^ — 1 arbitrary points 71, ... , 7jv-i according to interpolation formula 



Af-l N-l N-l 



m=l j=l j=l mj^j 

in our case 

N 



m = (a + A + E n - + E ^(^^o n (5-9) 



F(A) = n(A-A,.) (5.10) 

k=l 

with /i = — Ai — . . . — Ajv . Hence, we obtain ( p.la| ). 

Further, we consider the relation ( |3.1b| ). Using ( |2.13| ) one obtains 



„ N-l 

CUK)^Xu...,Xn = -e"^ / /""'(7)Q(7; A) JJ (A„ - 7,)^-r,Al+...+A^. (5.II) 

c ^=1 



Clearly, 

e"^^^,A,+...+A^ = ^^M+...+x^-in (5.12) 
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and, therefore, the relation ( 5.11|) acquires the form 



N-l 



C^.(A„)^Al,...,A^ = (-1)"^ [ /x-i(7)g(7; A) JJ(7, - \n)'^^M+-+XN~ihd'y (5.13) 

c 

Evidently, the function (5(7; A) satisfies to equation 

N-l 

n - >^n)Q{T, A) = 2^-'g(7; A - ^neO (5.14) 
Hence, we prove that function ( [4.3| ) obeys the relations ( |3.1b| ). 

The final step is to prove that the function ( |4.3| ) is the genuine Whittaker function. The inte- 
grand in ( [4. 3D decreases exponentially as ■jj —ioo , {j = 1, . . . , N—1) and, as consequence, 
the integrals over large semi-circles in the lower half-plane vanish. Using the Cauchy formula 
to calculate the integral (|473|) in the asymptotic region Xk+i ^ Xk, {k = 1, . . . , N — 1) , it is 
easy to see that the asymptotics of the function '0Ai,...,A]v determined precisely in terms 
of the corresponding Harish- Chandra functions (see, for example, [^): 

M^) = E /l-^-(-^'^)/^nr( '^\7^' )^^^^^'"^ + 0(max{e---^+^}^""j) (5.15) 

s&W j<k ~ 

(in the last formula the summation is performed over the Weyl group). Hence, we construct 
exactly the Weyl invariant Whittaker function. Moreover, using the Stirling formula for the 
r -functions, it is easy to see that the asymptotics ( p.2| ) hold. Theorem is proved. ■ 



6 The Mellin- Barnes representation 

Theorem 6.1 Let a set W'jjkW be the lower triangular N x N matrix. The solution to 
eqs. ( ^^) can be written (up to inessential numerical factor) in the form of multiple Mellin- 
Barnes integrals: 

^/'7^,,...,7^^(a^i,---,a;iv) = 

11 n TT 1 Xn[lnk " 7n-l,fc ) M 1 dljk 



n=l 



np / 7nj Ink \ 



n,k=l J j,k=l 



C 

j,k=l 
j<k 

(6.1) 

where the integral should be understand as follows: first we integrate on ju over the line 
Im7ii > max{Im72i, Im722} ; then we integrate on the set (721,722) over the lines Im72j > 
maXm{ Im73m} and so on. The last integrations should be performed on the set of variables 
(7iv-i,i • • • , 7iv-i,iv-i) over the lines Im7jv-i,fc > maXm{Im7jv,m} . 



The proof is straightforward resolution of the recurrent relations ( |4.3| ) starting with trivial 
Whittaker function 'ip^_^^{xi) = exp{-|7iiXi} . 
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